Pinning and depinning of a classic quasi-one-dimensional Wigner crystal in the 

presence of a constriction. 
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We studied the dynamics of a quasi-one-dimensional chain-like system of charged particles at low 
temperature, interacting through a screened Coulomb potential in the presence of a local constric- 
tion. Fhe response of the system when an external electric field is applied was investigated. We 
performed Langevin molecular dynamics simulations for different values of the driving force and for 
different temperatures. We found that the friction together with the constriction pins the particles 
up to a critical value of the driving force. The system can depin elastically or quasi- elastically 
depending on the strength of the constriction. The elastic (quasi-elastic) depinning is characterized 
by a critical exponent f3 ~ 0.66 (/3 ~ 0.95). The dc conductivity is zero in the pinned regime, it has 
non-ohmic characteristics after the activation of the motion and then it is constant. Furthermore, 
the dependence of the conductivity with temperature and strength of the constriction was investi- 
gated in detail. We found interesting differences between the single and the multi-chain regimes as 
the temperature is increased. 

PACS numbers: 71.10.-w, 52.65.Yy, 45.50.Jf 



I. INTRODUCTION 



In the last years the interest in mesoscopic systems 
consisting of interacting particles in low dimensions or 
confined geometries has seen a sustained growth. A class 
of quantum anisotropic systems exhibiting "stripe" be- 
havior appears in the quantum Hall regime |lj , in charge 
density waves (CDW)0|, in manganite oxides and high- 
T c superconductors [3j where strong electronic correla- 
tions are responsible for the formation of these inhomo- 
geneous phases. Another class of confined quasi one- 
dimensional (Q1D) geometries appears in different fields 
of research and some typical and important examples 
from the experimental point of view are: ordered elec- 
trons on microchannels filled by liquid Helium Q, mi- 
crofluidic devices 0, colloidal suspensions Q and con- 
fined dusty plasma |8(. On the atomic scale a chain- like 
system was found in compounds such as Hg 3 _ s AsFR 
and in low dimensional systems formed on surfaces [1C 
These kinds of interacting systems, which tend to form 
periodic or ordered structures when the density of par- 
ticles and the temperature are low enough (i.e. Wigner 
crystallization |lll 12]), can exhibit a remarkable vari- 
ety of complex phenomena when they are driven by an 
external force. Many of these phenomena, which arise 
from the interplay between periodicity, disorder, nonlin- 
earities and driving, are still poorly understood or even 
unexplored. For numerous such experimental systems, 
transport experiments 0, Q, 0] are a useful way to 
probe the physics (and sometimes the only way when 
direct methods, e.g. imaging, are not available). It is 
thus an interesting and challenging problem to obtain 
a quantitative description of their non-linear dynamics. 
One striking property exhibited by all these systems is 
pinning, i.e. at low temperature there is no macroscopic 
motion unless the applied force / reaches a threshold crit- 



ical value f c . There is a quite extensive literature about 
the dynamical properties near the depinning threshold 
El EE 111 , mostly in the context of CDW [Hill EH- 



The aim of this paper is to provide a description of 
the properties of a Q1D Wigner crystal in the presence 
of a local constriction potential, thermal noise and an 
external driving force. Most of the previous theoretical 
and experimental works are on moving 2D or 3D lattices 
and glasses (see Ref. [^ and references therein). The 
effect of confinement into a mesoscopic channel has not 
yet been deeply investigated. 

Our classical model is very ductile because of its scal- 
ability and it is suitable for the description of diverse 
confined systems, as electrons on liquid helium, colloids 
and complex plasmas. We should stress that the classical 
approach, which is naturally valid in the case of colloids 
and complex plasmas because of the microscopic size of 
the particles, is still valid for pure quantum objects as 
electrons when they exhibit Wigner crystallization. In 
the crystal phase the electrons become localized, and 
thus distinguishable. In this case the De Broglie thermal 
length \r> is much smaller than the interparticle spac- 
ing, hence the quantum aspect of the original fermionic 
problem does not play a crucial role and the classical 
treatment for the system is an accurate one. Several 
generic aspects of the present model without a constric- 
tion and in the linear regime were recently investigated 
by the present authors |23|. 

A narrow channel with a constriction can be readily 
realized in a colloidal system or in a dusty plasma. Ad- 
ditionally the problem we deal with could also be of 
interest in nanoscale wires. Flux-line-lattice flow has 
been studied in a novel superconducting devices contain- 
ing straight, nanometer-scale, weak-pinning channels in 
a strong-pinning environment j24[ . By introducing a con- 
striction into the weak-pinning channels many features of 
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the model we propose can be investigated experimentally. 

The paper is organized as follows: we first give in Sec. 
II an overview of the model and of the numerical meth- 
ods used. In Sec. Ill we describe the zero temperature 
phase diagram in the absence of any external driving 
force, stressing the differences in the ground state con- 
figurations near the constriction. Sec. IV is devoted to 
the study of the dynamics of the system, in particular 
we concentrate on the velocity vs applied driving force 
curves and on the conductivity. In Sec. V we discuss 
the interplay between driving force, thermal disorder and 
constriction potential, focusing on the difference between 
the single chain and multi-chain regime. In Sec. VI we 
comment on the analogies and differences with the case 
of moving lattices and glasses and other models in which 
quenched disorder or pinning potentials are present. Fi- 
nally we conclude in Sec. VII. 



II. MODEL AND METHODS 

The system consists of an infinite number of classical 
identical particles with charge q and mass m, moving in 
a plane with coordinates r— {x,y). The particles inter- 
act through a screened Coulomb (Yukawa-type) poten- 
tial, where the screening length A is an external param- 
eter. We impose a parabolic confining potential in one 
direction, namely in the y-direction, and a constriction 
potential with Lorentzian shape centered in the origin of 
the axes. The Hamiltonian of the system is given by: 
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independent of the specifics of the system under investi- 
gation. The dimensionless time and temperature, which 
are essential quantities in what follows, are respectively 
T' = T/[k B (muj 2 q 4 /2e 2 ) 1 / 3 } and r = tcj . 

The zero temperature configurations for different den- 
sities, namely different number of particles, and different 
values of the parameters in the Hamiltonian were calcu- 
lated by the Monte Carlo (MC) technique using the stan- 
dard Metropolis algorithm as it was done in Ref. |2^| . We 
have allowed the system to approach its equilibrium state 
at some temperature T, after executing 10 5 -j- 10 6 Monte 
Carlo steps. In order to reach the T = equilibrium 
configuration the technique of simulated annealing was 
used: first the system was heated up and then cooled 
down to a very low temperature. We introduced peri- 
odic boundary conditions along the x-direction in order 
to simulate an infinite long wire. Typically a simulation 
cell of length L = 100 (in dimensionless units) centered 
around the origin of the axis was used. This choice was 
motivated by the fact that for larger L the changes in 
static and dynamical properties are negligible, especially 
for a' > 0.5 and k > 1. 

Because of the presence of the constriction we found 
many more metastable states, which complicates the nu- 
merical approach. We will elaborate on this point in the 
next section. 

After reaching the T — equilibrium configuration, we 
introduced an external electrical field in the x-direction, 
or in other words we considered the effect of an external 
driving force / and calculate the transport properties of 
the system. We also considered the effect of temperature 
and thermal noise, c oup ling the system to additional de- 
grees of freedom [25|, |2g or to a heat bath. The Langevin 
dynamics |27| is the most appropriate one to include such 
effects. The Langevin equations for the x and y compo- 
nents of motion in dimensionless units are respectively: 



where e is the dielectric constant of the medium the par- 
ticles are moving in, loq measures the strength of the con- 
fining potential, Vq is the maximum of the potential of 
the constriction which has a full width at half maximum 
2/ a. Introducing a suitable system of units, the Hamil- 
tonian can be rewritten in a dimensionless form. Wc 
define: r = (2q 2 /meLo 2 ) 1 / 3 and E = {muj 2 q A /2s 2 ) 1 / 3 
as unity of length and energy, respectively. After using 
dimensionless units the Hamiltonian takes the form: 



_ exp(-K|r^ -r'j\ 
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where H' = H/E Q , k = r /A, ~r' = ~?/r , V ' = V /E 
and a' = r^a. 

This transformation is particularly interesting because 
now the Hamiltonian no longer depends on the mass of 
the particles, the dielectric constant of the medium and 
the frequency of the parabolic confinement, that is it is 
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where 7 is the friction coefficient, which is an external pa- 
rameter as well as k, and £(T') is a random force, repro- 
ducing the thermal noise, with zero average and standard 
deviation 



(^ir(T)^s(cr)) = ^T'5 l j5 rs 8{T - a). 

where (r, s) = (x, y). The driving force / and the random 
force £ are measured in units of mw^Q, while the fric- 
tion coefficient 7 is measured in units of ojq. We used the 
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same simulation cell and the same boundary conditions 
as in the case of the MC simulations. It should be no- 
ticed that in the case of colloids or vortices in type II 
superconductors |29| the motion is overdamped and Eqs. 
(3a) and (3b) are simplified: the hydrodynamic interac- 
tions can be neglected and Eqs. (3a) and (3b) reduce to 
a system of coupled first order differential equations. 

We considered here the more general problem, includ- 
ing also the hydrodynamic terms. In order to integrate 
the equations of motion we used a quasi-symplectic algo- 
rithms of "leap frog" type [3(J in the form: 



. . At . . 
n = n{t) + —Vi(t), 



Vi(t + At) 



C-2 



c 1 v l (t) + At—(fi) + dir/t 
At 



r t (t + At) = f t + — Vi{t + At), 

where At is the time step, and rji are Gaussian variables 
with standard deviation equal to 1 and average equal to 
0; the constants c\, C2 and d\ are given respectively by: 
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It was shown that this integration scheme has a good 
stability and runs rather fast, furthermore it is well be- 
haved in the limit 7 — > [3lJ • 

The driving force was increased from zero by small 
increments. A time integration step of At = 0.001 was 
used and averages were evaluated during 2 x 10 5 steps 
after 2 x 10 6 steps for equilibration. 



III. GROUND STATE CONFIGURATIONS 

The ground state configuration is the result of compet- 
itive effects, that is the electrostatic repulsion, the con- 
fining potential that tries to keep the particles close to 
the x axis and the Lorentzian constriction potential that 
prevents the particles from settling close to the y axis. In 
Fig. 1 the contour plots of the sum of both potentials for 
two different values of V ' are shown. Depending on the 
values of (increasing) Vq and (decreasing) a' the saddle 
point at (x,y) = (0,0) becomes more pronounced. 

As discussed in Ref. 23], in the absence of any con- 
striction the charged particles crystallize in a number of 



chains. Each chain has the same density resulting in a to- 
tal one-dimensional density n e . If a is the separation be- 
tween two adjacent particles in the same chain, it is pos- 
sible to define a dimcnsionlcss linear density n e = lro/a, 
where I is the number of chains. In the case of multi- 



ple chains, in order to have a better packing (or in other 




FIG. 1: Contour plots of the sum of the confinement and con- 
striction potentials for two different values of the parameters: 
(a) Vi = l,a' = 1, (b) V ' = 5,a' = l. 



words to minimize the interaction energy by maximizing 
the separation among particles in different chains), the 
chains are staggered with respect to each other by a/2 in 
the x-direction. For low densities the particles crystal- 
lize in a single chain; with increasing density a "zig-zag" 
(continuous) transition [32l | occurs and the single chain 
splits into two chains. Further increasing the density the 
four-chain structure is stabilized before the three-chain 
one. This first four-chain configuration has a relatively 
small stability range after which it transits to a three 
chain configuration. For higher values of the density, the 
four-chain configuration attains again the lowest energy. 
Then a further increase of n e will lead to more chains, 
that is five, six and so on. The structural transitions 
are all discontinuous (i.e. first order), except the 1 — * 2 
transition. 
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In the presence of the constriction potential the ground 



state configurations are modified near the constriction 



4 




20 16 -C S 6 10 - 5 20 25 



09 

0.5 
04 
02 
>. 03 
•02 
■'j i 



-0.3 







H 










i 




' N=ieo 








* 




^ 



-» -ao .is .10 •& o 5 id 15 -jo 



1 2 
03 
04 
> j 
-0 4 
-03 



•*«••••• 



fe 



,* » . . 



,*■■*»*■* 



' 2 25 



U 20 26 



I 

Od 

0.4 
>> 0.0 
•04 

■ 3 

'2, 

20 
IS 
10 
05 
>> 00 
-0 5 
-10 
-15 
-20 











(s) 


■ N»240 














■ ■ 

t • 
















25 -20 -15 


■l" 


-5 C 5 




15 20 25 










0) 






• 








mm 


''•Mill 


•MM 






IHtf 


• 

!•*».,* Illllll 














u=1 



-25 -20 -15 -1C -5 C 5 -0 15 



04 

o: 

-02 
•04 



(b) 



! 



u'=01 



:c t6 10 5 



'C II 20 25 




•30 OS .10 ■) i K 15 a 8 



I 



■ - 

C4 

CO 
-C4 





-04 
; 3 U.*t*. 



N=22« 



(i-0 1 

,*»»•■•»»•••••*•*♦*•»» 



•••••I 



••*•. 



I5 



15 2C 25 



8 
- . 



-25 -20 -15 -10 -5 5 13 1= 2C 25 




•• ii 



(1=0 1 



•25 -20 -16 -10 -5 5 10 16 20 25 
X 



FIG. 2: The ground state configuration in the center of the simulation cell for k = 1, Vq = 1 and different values of the number 
of particles TV in a simulation cell of length L — 100 and two different values of the constriction width: (a) TV = 80, a' = 1, (b) 
TV = 80, a' = 0.1, (c) TV = 160, a' = 1, (d) TV = 160, a' = 0.1, (e) TV = 228, a' = 1, (f) TV = 228, a' = 0.1, (g) TV = 240, a' = 1, 
(h) TV = 240, a' = 0.1, (i) TV = 440, a' = 1, (j) TV = 440, a' = 0.1. In the case of a' = 1 the effect of the constriction potential 
is significant only around a narrow region around (x,y) = (0,0), while for a' = 0.1 it is appreciable all along the lenght of the 
simulation cell. It is interesting to notice that for a' = 1 the reentrant behavior of Ref. |23j|. with the four chain arrangement 
stabilized before the three chain arrangement, is still present, while it is absent for a' — 0.1. For small values of a' the 
configurations are highly inhomogeneous. 
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(see Fig. 2), but the particles are still organized in chains 
far away from this constriction. Close to the saddle point 
the particles do not arrange themselves in ordered chains. 
The particles near the constriction lead to a significant in- 
crease of the number of metastable states. Consequently 
the procedure of simulated annealing has to be more ac- 
curate than in the case of the absence of a constriction, 
which means that several intermediate temperature steps 
have to be considered. Sometimes the MC simulations do 
not provide us with the "exact" ground state, as it is seen 
for instance in Fig. 2 (c), (f) and (j), where the final con- 
figurations are not perfectly symmetric with respect to x 
and y, while this is expected because of the symmetry of 
the Hamiltonian. 

When the full width at half maximum of the constric- 
tion potential is short enough (a' > 0.5), or in other 
words the effect of the constriction is significative only in 
a narrow region around x = 0, it is still possible to de- 
fine a local density because the system exhibits a homo- 
geneous spacing among charged particles except in the 
vicinity of the saddle point. Thus, excluding these re- 
gions, it is still meaningful to consider n e — lr$/a, where 
I is the number of chains. In this case the same chain ar- 
rangements, i.e. l-^2^4^3^4— >5 and so on, as 
in the case where the constriction is absent, is found (see 
Fig. 2 (a), (c), (e), (g) and (i)) with increasing density, 
but with the difference that all the structural transitions 
are now discontinuous. 

For smaller values of a, that is for larger interaction 
ranges of the constriction potential, the system is highly 
inhomogeneous and even shows coexistence of different 
chain phases (see Fig. 2 (d), (f), (h) and (j)). In a 
certain sense the constriction introduces a local disorder 
into the system. 

In order to make these affirmations clearer, we plot in 
Fig. 3 the inverse interparticle spacing, i.e. the density, 
in the single chain configuration as a function of the dis- 
tance from the origin of the coordinates for different val- 
ues of k and a' and Vq = 1. It is evident that the density 
is an increasing function of the inverse screening length 
k, because the particles can stay closer together as the 
electrostatic repulsion is weaker. For small values of a' 
the density is an increasing function of the distance over a 
large range of a;- values (Fig. 2(a)), while for a' — 1 (Fig. 
2(b)) this range can become very small and the density 
becomes very quickly independent of x. Notice that for 
x — > oo the chain density should become independent of 
the parameters of the constriction. 



IV. DYNAMICAL PROPERTIES 



ra 



0.80 
0.75 
0.70 
0.65 
0.60 
0.55 
0.50 
0.45 

0.75 
0.70 
0.65 
0.60 
0.55 
0.50 



a ' =0 - 4 /' K =5 

N=80 / K=1.0 

K=0.1 



1 chain structure 



(a) 



10 



20 30 
X 



40 



50 



V'=1 


<x'=1 y''' 




— K=5 


. N=80 


— K=1.0 


* 

' ...... 

1 
t 


•~K=0.1 


1 chain structure 
i.i.i. 


(b) 

■ 



10 



20 



30 



40 



50 



FIG. 3: The density as a function of the distance from the 
origin in the single chain configuration for 80 particles, for a 
large (a) and short (b) interaction range of the constriction 
potential. Because of the symmetry of the system around 
y — only the results in the right part of the simulation cell 
are reported. 



tion respectively, encountered in complex plasmas |35|| . 
We also fixed the value of a' = I, that means that we 
deal with short range constrictions. 

The first obviously important quantity to determine is 
the velocity v' as a function of the applied force /. In 
the absence of thermal fluctuations, i.e. T — 0, and in 
the absence of the constriction potential, Eqs. (3a) and 
(3b) become: 



dr 



-y— 

2 ^ dx' 

3 1 



d exp(-K|r^ - r'A) 



+ /; (4a) 



When a constant electrical field E is applied to the sys- 
tem in the ^-direction, it produces a longitudinal driving 
force /. The charged particles then are pushed along the 
direction of the driving force. In what follows we consider 
mainly systems for which k = I and 7 = 0.2 , which are 
typical values for the inverse screening length and fric- 



dr 
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Furthermore, because in the equilibrium configuration 
the net force acting on every particle, due to electrostatic 
repulsion and confinement, is zero, that is 
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Eqs. (3a) and (3b) can be ulteriorly simplified and one 
obtains the uncoupled equations: 
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whose stationary solutions are respectively: 
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This shows that in the absence of thermal noise and 
constriction the total effect of the external driving force 
is a sliding of the ordered structure with a drift velocity 
which is directly proportional to the driving force and 
inversely proportional to the friction. More in general, 
when the leading term in the equation of motion is the 
driving force, one should expect that the drift velocity is 
v' x — f /"/, or in other words that the system behaves like 
a classical two-dimensional Drude conductor (3^1 . This 
feature has been observed in experiments Q and in nu- 
merical simulations |34| . In the presence of a constriction 
and thermal noise, v' x is no longer a linear function of the 
(6a) driving force, as we will discuss in the next subsections. 



(5a) 



(5b) 



A. Pinning 

The system is pinned until the applied driving force 
reaches a threshold value f c . The pinned structures in 
the presence of the driving force show substantial differ- 
ences with the ground state configurations in the absence 
of a driving force. The particles move into the direction of 
the driving force and accumulate in front of the constric- 
tion, that exerts a force, which is opposite to the driving 
force. If f < f c , new static configurations are reached, in 
which the electrostatic repulsion and the repulsive force 
due to the constriction balance the driving force. The 
situation is depicted in Fig. 4, where the driving force is 
in the positive direction of the x axis. In the case of low 
constriction barrier height, the chain structures are rela- 
tively homogeneous, although the inter-chain distance is 
smaller at the left than at the right of the constriction 
because of the external drive. In the case of high con- 
striction barrier, the chain structures are no longer ho- 
mogeneous. As can be seen in Fig. 4, different number 
of chains can coexist in the same configuration. Since 
the energy to overwhelm the barrier is quite large, the 
particles tend mainly to accumulate at x < 0, which pro- 
duces a density gradient and consequently a splitting into 
a larger number of chains where the density is larger, be- 
cause obviously in such a case the electrostatic repulsion 
among particles is larger. 

It should be noticed that the nature of pinning for the 
system under investigation is different from the pinning 
often studied in literature for e.g. colloidal systems and 
vortex lattices. In these cases the pinning is the result 



of some kind of disorder (in most cases quenched disor- 
der) or, in other words, the effect of the substrate. It 
is introduced into the system and modeled as randomly 
placed point-like pinnin g ce nters producing an attractive 
gaussian potential |2J, 1361 |3jJ |3J1 |3j| or as randomly 
placed parabolic traps 28]. In our case, the pinning is 
the effect of a constraint, the particles have not enough 
energy to overcome the constriction barrier and conse- 
quently there is no net motion. We also investigated the 
case of negative Vq, i.e. a single Lorentzian potential 
well. In that case we still observed pinning, but without 
the formation of highly inhomogeneous structures or, in 
other words, without the accumulation of many particles 
in the direction of the driving force. In the case of a neg- 
ative Lorentzian potential the same chain configurations 
are preserved along the simulation cell length even in the 
case of a large depth of the potential well. In the case 
of CDW the pinning potential can be either attractive 
or repulsive, indeed the the sign of the potential can be 
converted by changing the phase of the CDW by tt. In 
what follows we will limit ourselves to positive Vq. 

In Fig. 4, the trajectories of the particles are reported 
for a temperature T' = 0.002, well below the melting 
temperature. It is interesting to study for a fixed num- 
ber of particles and for a fixed temperature how by in- 
creasing the driving force the configurations change. The 
variation of the density along the constriction is shown 
in Fig. 5 for different values of the external driving force 
for a constriction height of Vq = 5 and width a' = 1. 
Increasing the driving force more and more particles ac- 
cumulate to the left of the constriction barrier in the di- 
rection of the driving force, corresponding to larger and 
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FIG. 4: Typical trajectories of the particles when the system is pinned for different number of particles N in a simulation cell 
of length L = 100 and different values of the height of the constriction barrier. The plots are for a temperature T' = 0.002, 
k = 1 and a' — 1. 
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FIG. 5: The density as a function of the distance along the 
simulation cell for different values of the driving force and for 
a temperature T = 0.002. 



larger densities. The density n e has a discontinuity at the 
constriction. For low values of the driving force, except 

I 



The first interesting observation is that the driven sys- 
tem does not break up into pinned and flowing regions, 
as observed in experiments and simulations of supercon- 
ducting vortices [4pL [III 142^ or colloids [2£|, or, in other 
words, the chain-like system under investigation does not 
exhibit plastic depinning. Once the driving force over- 
whelms the critical threshold / c , all the particles move 
together. The depinning can be either elastic or quasi- 
elastic depending on the height of the constriction bar- 
rier. In the case of a low barrier (V ' — 1 in Fig. 6) 
the particles move such that they keep the same neigh- 
bors, thus the system depins elastically. In contrast in 
the case of a high barrier (Vq = 5 in Fig. 6) a complex 
net of conducting channels is activated and the particles 
move without keeping their neighbors, that is the de- 
pinning is quasi-elastic. The quasi-elastic depinning is a 
feature closely related to the low dimensionality of the 
system and to the fact that we are considering a single 
constriction. We found that quasi-elastic depinning ap- 



in the vicinity of the constriction, n e is almost constant. 
But for larger values of the driving force it is always an 
increasing function of the distance along the simulation 
cell, except close to the constriction. 

Increasing T", we observed larger and larger oscilla- 
tions of the particles until the system is melted. As 
already reported in Ref. |2^|, also in this case we ob- 
served larger oscillations in the x-direction than in the 
y-direction, which is evidently a combined effect of the 
confining potential and the nature of the interparticle in- 
teraction. In the case of high temperature, close to the 
solid-liquid transition, and mainly in the case of large 
Vq the arrangements of the particles are slightly different 
from the ones shown in Fig. 4. Because high values of 
Vq in combination with the driving force produce a den- 
sity gradient in the chain-like structures, the melting is 
not homogeneous, with the coexistence of solid and liq- 
uid regions. It is beyond the aim of the present paper to 
discuss how the driving force induces local melting of the 
system. 

The critical force f c is evidently a function of the tem- 
perature T", it decreases with increasing temperature, 
that is the thermal motion aids the net motion of the 
particles. The critical force is also a function of the den- 
sity, i.e. the number of particles. In our simulations we 
observed that for larger densities f c becomes smaller. 



B. Depinning 

When the driving force / is larger than the threshold 
f c , the system exhibits Q1D flow. In Fig. 6 some typical 
trajectories of the depinned particles are reported, for 
T" = 0.002 and different values of / just above f c , 

I 

pears when the strength of the constriction potential is 
increased. In other infinite 2D systems of driven particles 
or vortices, where several pinning centers are present, a 
crossover from elastic to plastic depinning with increas- 
ing strength of the pinning potential occurs (see e.g. Ref. 
[28|). We will focus on this problem in the next subsec- 
tion. 

The region after the constriction barrier, as it is ev- 
ident from the trajectory patterns, shows features that 
deserve a deeper investigation in the case of high values 
of Vq. In Fig. 6 for V ' = 5 we found that at the right 
of the constriction some noise is present and the parti- 
cles flow disorderly. In order to explain this behavior we 
investigated the distributions of the x and y components 
of the velocity in narrow strips along the simulation cell 
length. We concentrated our attention on a system of 
N = 400 particles at V = 0.002, with k = 1, V ' = 5 and 
a' = 1. The results are reported in Fig. 7. 
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FIG. 6: Typical trajectories of the particles when the system is depinned for different densities and different values of the height 
of the constriction barrier. The plots are for a temperature T' = 0.002, k = 1 and a' = 1. 
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It is evident that v y is always normally distributed 
with average equal to zero, as expected, because it re- 
ceives contributions mainly from thermal noise, which is 
gaussian. The distribution of v x is still gaussian, but cen- 
tered around a value < v x >^ because of the external 
driving force, except in the neighborhood of the constric- 
tion where the strong interaction with the barrier gives 
a non-gaussian profile to it. What is interesting is the 
fact that the velocity above the depinning threshold has 
a pronounced gradient in the x-direction. From Fig. 7, 
it is clear that approaching the constriction from the left 
side the particles are slowed down; they receive a sudden 
acceleration when they pass the constriction barrier, then 
the velocity has a maximum in the right neighborhood 
of the constriction and finally it slows down again when 
approaching the edge of the simulation cell. In order to 
explain these highly non linear features, it is helpful to 
look at the profile of the force due to the constriction 
potential (see Fig. 8). This force has a significative mag- 
nitude only in a narrow region around x = 0. For x < 
it acts oppositely to the driving force, while for x > 
it enhances the driving force. There are two maxima for 
the intensity of the force located at x = ±l/y/3 and it 
is zero at the origin of the axis. Therefore, when the 
particles approach the constriction they start to feel this 
decelerating force and slow down. Because the system 
is strongly interacting the deceleration is seen not only 
in the left neighborhood of the barrier, but in a wider 
region. At x = the force is zero, for x > close to 
the constriction the force quickly increases and adds to 
the driving force, so the particles experience a sudden 
acceleration which produces a large velocity. After that 
the particles are accelerated only by the driving force and 



start to feel the effect of the particles on the opposite side 
of the simulation cell because of the periodic boundary 
conditions, so the velocity decreases again. In the case 
of low constriction barriers or large driving forces the x 
components of the velocity are more homogenously dis- 
tributed along the simulation cell. From the width of the 
velocity distribution it is possible to define an " effective 
temperature". According to the equilibrium probability 
factor 



P oc exp 



rnvz 



2k B T; 



+ 



eff 



2k B T* 



eff 



where T^jj and take into account also the contribu- 
tions due to the potential energy, the fluctuations of the 
velocity components are related to the effective temper- 
ature by: 



T* ff = m<(v' x -<v x >f >/k B , 
T y eff = m< (v' y - < v' y >) 2 >/k B . 

In our dimensionless units and our specific case this 
yields: 



1 eff 



2 < (v' x - < v' x >) 2 >, 
T'y f = 2<v' y 2 >, 



respectively. The calculated effective temperatures are 
reported in Table I. 



It is worth to notice that the effective temperatures are 
the same as the simulation temperature in the regions 
far away from the constriction barrier where the velocity 
fluctuations are nicely described by a normal distribu- 
tion (see Fig. 7). The effective temperatures and 
7^// increase when approaching the constriction. In the 
strips -0.5 < x < 0.5 and 1.5 < x < 2.5 both T'* f} and 

T'gff T' . This is evidently a result of the strong inter- 
action with the barrier which increases significantly the 
fluctuations in the velocity. Actually, the spreading of 
the distribution of the velocities is one to two orders of 
magnitude larger in the constriction region with respect 
to the regions where the constriction potential is almost 
zero. 





rplx 
1 eff 


1 eff 


-40.5 < x < -39.5 


0.0018 


0.0020 


-20.5 < x < -19.5 


0.0020 


0.0021 


-2.5 < x < -1.5 


0.0060 


0.0051 


-0.5 < x < 0.5 


0.28 


0.031 


1.5 < x < 2.5 


0.14 


0.26 


19.5 < x < 20.5 


0.0024 


0.0023 


39.5 < x < 40.5 


0.0019 


0.0021 



TABLE I: Effective electron temperatures corresponding to 
the situation of Fig. 6(h) in the different regions studied in 
Fig. 7. 



In the elastic regime the velocity fluctuations could be 
fitted to a Gaussian distribution both for v x and v y , even 
in the vicinity of the constriction. Around the barrier, 
the strong interaction effect is felt as an increase of the ef- 
fective temperature, which is much less pronounced than 
in the quasi-elastic regime. 



Regarding to the noise observed in the trajectory plots, 
it is essentially due to the fact that the particles merge 
from the constriction with a relatively large y component 
of the velocity. As is seen form Fig. 7(d), the distribu- 
tion of the v y spreads over a quite large range of values 
while passing the constriction. Indeed, the standard de- 
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FIG. 7: The v x and v y distributions along the simulation cell length for a system of N = 400 particles, for n = 1, Vo = 5 and 
a' = 1 and for a temperature T' = 0.002 and driving force / = 0.05. The superimposed curves are the normal distribution 
curves generated using the mean and standard deviation of the data. 



12 




, i , i i i , i , i , i i i , I ■* 

■*00C\I^-Ot-CMC5^- 
i i i i 



80J0; UOjpUJSUOO 



FIG. 8: The constriction force profile as a function of x for 
Vq = 5 and a' = 1 



viation of the v y distribution around x = is one order 
of magnitude larger than in the other regions, as already 
mentioned. This feature is a consequence of the fact that 
very close to the barrier the constriction force is very 
small (it is zero at x — 0), the particles proceed slowly 
and, thereby, they strongly undergo the effect of the con- 
fining potential which accelerates them in the y-direction, 
producing a significant y component of the velocity. This 
is also confirmed by the fact that passing the constric- 
tion barrier a narrowing of the chain structures is always 
present (see Fig. 6). Notice from Fig. 7(d,e) that the 
velocity fluctuations are no longer described by a nor- 
mal thermodynamic equilibrium distribution and that in 
particular for v x there are large deviations. 



As first predicted by Fisher the elastic depinning ex- 
hibits criticality an d the velocity vs force curves scale 
as v — (/ — f c )P. This scaling has been extensively stud- 
ied in 2D CDW systems where f3 = 2/3 @ E^. It is, 
however, still an open issue whether this exponent is the 
signature of an universality class and whether it depends 
on the particle-particle interaction. Actually, in other 
investigations on elastic depinning of driven colloidal lat- 
tices the findings were (3 ~ 0.5 ||g and (3 = 0.92 ± 0.01 
|34j . Other studies on plastic depinning with filamentary 
or river-like flow have shown a velocity-force curve scal- 
ing with ^ = 2.2 0J = 1.94 ± 0.03 for colloids, 
(3 — 2.0 for electron flow simulations in metallic dots |47j 
and (3 = 2.22 for vortex flow in superconductors |3^ |. 

As pointed out by Le Doussal and Giamarchi, for an 
infinite size 2D system true elastic depinning is not ex- 
pected since dislocations and defects, acting as pinning 
centers, should appear at large scales [22J. Both the sim- 
ulations and the experiments are, however, always for 
finite size systems and consequently elastic depinning is 
possible where the distance between dislocations may be 
larger than the system size. 

In Fig. 9, we report the v — f curve in the case of elas- 
tic and quasi-elastic depinning for different number of 
particles, i.e. for different chain arrangements. It should 
be noticed that the critical exponent does not depend on 
the number of chains in the system. For all the investi- 
gated chain configurations we obtained on average that 
(3 ~ 0.66 in the case of homogeneous channel flow, that is 
elastic depinning, and (3 ~ 0.95 in the case of inhomoge- 
neous channel flow, that is quasi-elastic depinning. For 
negative value of Vq we found similar values for the criti- 
cal exponent in the elastic and quasi-elastic regime. The 
value of the critical exponent could, therefore, be consid- 
ered as a clear signature of the kind of depinning. Our 
results are consistent with the findings on CDW systems 
and colloids, mentioned before. With increasing temper- 
ature but below the melting temperature, we observed a 
broadening of the conducting channel or some changes in 
the structure with some chains collapsing (we will pro- 
vide more details about this point in the next subsection), 
but no significant dependence of the critical exponent on 
temperature was found (within our fitting errors). 



The question of whether in confined systems there is an 
universal exponent for elastic and quasi-elastic depinning 
cannot be answered conclusively. We found that the crit- 
ical exponent are not affected by the value of k, as it can 
be seen in Table II: going from k = 0.2 (nearly Coulomb 
interaction) up to k = 5 (short range interaction), the 
value for the critical exponent stays the same. Although 
we found that the critical behavior is independent of the 
particle-particle interaction, further investigations with 



other kind of pinning potentials are required in order to 
affirm that the elastic, quasi-elastic and plastic depinning 
belong to different universality classes. 

In our simulations neither in case of elastic nor quasi- 
elastic depinning history dependence was found. The ve- 
locity vs applied drive are not hysteretic, and we obtained 
the same result for increasing and decreasing values of /. 

Finally, in the case of very small a, that is for a wide 
constriction interaction, because of the density gradient 
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FIG. 9: Velocity v x vs applied drive / for the elastic depinning ((a) and (b)), for the quasi-elastic depinning ((c) and (d)) for 
K = 1 and a' = 1. The dashed curves are the best fitted power law behavior for T' — 0. At zero temperature the curves 
exhibit a sudden jump at the depinning threshold, while at finite temperature they are smoother. The critical exponents are 
independent of the density, i.e. the number of chains. 



producing the coexistence of different chain structures, 
we always observed quasi-elastic depinning. 



C. Crossover from elastic to quasi-elastic depinning 

It is interesting to investigate the values of (3, or in 
other word the kind of flow, as a function of Vq. For 
colloidal systems a sharp crossover from elastic to plastic 
depinning was found with increasing strength of the sub- 
strate disorder, accompanied by a sharp increase in the 
depinning critical force [28|. Carpentier and Le Doussal 
studied theoretically the effect of quenched disorder on 
the order and melting of 2D lattices and found a sharp 
crossover from the ordered Bragg glass (where there are 
no defects) to a disordered state |43| • They also predicted 
that the depinning threshold increases at the order to dis- 
order transition due to the softening of the lattice, which 
allows the particles to better adjust to the substrate. A 
similar mechanism could account for the peak effect ob- 
served in vortex matter in superconductors |49| , in which 
the depinning threshold rises abruptly with increasing 
applied magnetic field. 

We found a crossover from elastic to quasi-elastic de- 



pinning as the barrier height of the constriction is in- 
creased. This is analogous to the crossover from the elas- 
tic to plastic flow encountered in other systems. As can 
be seen in Fig. 10, the behavior of j3 as a function of 
Vq is almost step-like, and the crossover takes place in a 
narrow range of V ' values. We also observed increasing 
values of the critical threshold f c . It is beyond the scope 
of this paper to determine whether the elastic to quasi- 
elastic crossover is a first or second order transition and 
how temperature influences this transition. However, the 
relative smoothness of the curves in Fig. 10 suggests a 
possible second order transition. Furthermore, the effect 
of increasing temperature should reasonably result in a 
shift to lower values of f c or V Q ' for the transition from 
elastic to quasi-elastic flow. 

It is evident from Fig. 10 that the crossover shifts to- 
wards lower values of the potential barrier as the inverse 
screening length is increased, which means that the par- 
ticles with stronger interparticle interactions can flow in 
a more ordered way. 

The values of the critical exponent (3 shown in Fig. 10 
can be fitted by the curve atanh (b(V ' — V Q )) + c, where 
a, b, c and V are the fitting parameters. In particu- 
lar, the sum a + c gives the value of (3 for the case of 
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elastic 
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quasi-elastic 


2 


0.25 


0.65 ±0.04 


elastic 


2 


1 


0.66 ±0.03 


elastic 


2 


5 


0.96 ±0.03 


quasi-elastic 


3 


0.25 


0.65 ±0.03 


elastic 


3 


1 


0.95 ±0.06 


quasi-clastic 


3 


5 


0.97 ±0.05 


quasi-elastic 


4 


0.25 


0.68 ±0.04 


elastic 


4 


1 


0.94 ±0.04 


quasi-elastic 


4 


5 


0.91 ±0.05 


quasi-elastic 


5 


0.25 


0.67 ±0.04 


elastic 


5 


1 


0.98 ±0.06 


quasi-elastic 


5 


5 


1.02 ±0.08 


quasi-elastic 



TABLE II: The critical exponent (3 for different values of the 
inverse screening length k and the constriction barrier height 



K 


a 


b 


c 


% 


0.2 


0.144 ±0.007 


2.1 ±0.6 


0.799 ± 0.007 


3.28 ± 0.06 


1 


0.131 ±0.005 


2.4 ±0.4 


0.790 ± 0.004 


2.90 ±0.05 


1.5 


0.139 ± 0.005 


3.2 ± 0.4 


0.796 ± 0.005 


2.12 ±0.06 


2 


0.143 ± 0.006 


4.2 ± 1.2 


0.813 ± 0.006 


1.78 ±0.04 


3 


0.152 ±0.005 


4.8 ± 1.3 


0.811 ±0.007 


1.24 ±0.05 


4 


0.146 ±0.007 


5.2 ±1.3 


0.802 ± 0.008 


0.94 ± 0.06 


5 


0.160 ±0.008 


5.4 ± 1.4 


0.80 ±0.01 


0.61 ±0.03 



TABLE III: The fitting parameters for the crossover curves 
for different values of the inverse screening length k. 

quasi-elastic depinning, while the difference c — a gives 
the value of j3 for the case of elastic depinnig, V can 
be identified with the value of the constriction height for 
which the crossover from elastic to quasi-elastic depin- 
ning takes place and b is related to the sharpness of the 
transition or, to be more precise, it is related to the in- 
verse width of the transition region. The results of the 
fits are reported in Table III. 

As it can be seen from Table III, the constriction height 
for which the quasi-elastic regime is established is a de- 
creasing function of the inverse screening length k, while 
the sharpness of the transition is an increasing function 
of k (see also Fig. 11). Again it confirms that for long 
range interactions the particles can flow more orderly. 
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FIG. 10: The critical exponent /3 as a function of the con- 
striction barrier height Vq for different values of the inverse 
screening length: (a)« = 0.2, (b) k = 1 and (c) k — 5. With 
increasing values of n the crossover region shifts to lower val- 
ues of Vq. The solid lines are guides to the eye. 

The errors in the fitting parameters are small for a, c 
and V and relatively large for b. 

In Fig. 11, the values of V and b as a function of k 

are reported. The curve V = V q (k) is well fitted by 
a Lorentzian l/(p ± qn 2 ) with p = 0.302 ± 0.005 and 
q = 0.058 ± 0.002 and the inverse width (see inset of Fig. 
11) by the Pade approximation (g + hn 2 )/(l + ln 2 ) with 
g = 1.9 ± 0.2, h = lA± 0.4 and I = 0.23 ± 0.07. 

We should stress that the physics behind the crossover 
from clastic to quasi-elastic depinning is different from 
the case of quenched disorder, where with increasing dis- 
order strength the ordered structure is softened and par- 
ticles can better adjust to the substrate. In our case 
the accumulation of particles in the vicinity of the con- 
striction barrier and their mutual repulsion give rise for 
high values of Vq to a complex arrangement of conducting 



FIG. 11: The transition points Vq as a function of the inve: 
screening length k. The dashed line is the fitting curve l/(j 
qn 2 ). The inset shows the behavior of the inverse width 
the transition region as a function of k. The dotted line is i 
fitting curve (g + hn 2 )/(l + Ik 2 ). 



channels, in which the nearest neighbors of each particle 
change, i.e to the impossibility of elastic flow. 



D. Conductivity 

According to Ohm's law the current density of a clas- 
sical system of charged particles is proportional to the 
applied electric field: 



= a-E, 



3=o 



(7) 



where a is the specific conductivity. The conductivity can 
be in general expressed as a second rank tensor. Because 
of the geometry of the investigated system and because 
the driving is in the x-direction, we are interested only 
in a xx , which we will refer in what follows as simply a. 

From the definition of j and from Eq. (6a) it follows, 
in dimcnsionlcss units: 



1 > 



(8) 



where a' = a/ (q 2 /mtuoro). As the definition of total 
density is not always an accurate one for our system, as 



FIG. 12: The ratio between velocity and driving force as a 
function of the drive strength for different values of the friction 
coefficient. The conductivity is non-Ohmic in a narrow region 
above the depinning threshold. The inset shows that for larger 
values of / Ohm's law is fulfilled with the conductivity as a 
constant. The curves are guides to the eye. 



mentioned above, we investigated the ratio v' x / f, which 
is directly proportional to a' through n e and which is a 
constant equal to I/7', according to Eq. (8). In Fig. 12 
the results of our calculations are reported for different 
values of the friction. 

When the particles are pinned the conductivity is ob- 
viously zero. It is interesting to notice that after the 
depinning threshold, there is a narrow region where the 
conductivity shows non-Ohmic features, going from zero 
to the saturation value I/7'. Afterwards when the drive 
is the leading effect in the equations of motion, the par- 
ticles behave as a classical Ohmic conductor. This be- 
havior is independent of the number of particles and of 
the height of the constriction barrier. For higher values 
of Vq the non-Ohmic conductivity region is enlarged. 

Studying the conductivity as a function of temperature 
some interesting features were observed. We investigated 
a rather wide range of temperatures from 0.001 to 0.018. 
From Ref. (2^ we know that for the single chain con- 
figuration in the absence of driving and constriction po- 
tential, the melting temperature is arbitrarily low, while 
in the multiF-chain configuration it is T' m ~ 0.015. The 
results of our calculations, for a driving force / = 0.05 
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FIG. 13: Dependence of the conductivity on the temperature, 
for (a) a single-chain structure and (b) a multi-chain struc- 
ture. For weak values of the potential barrier the conductivity 
of the single-chain structure is an increasing function of tem- 
perature, in all the other cases it shows a decreasing trend. 



FIG. 14: The average configurational energy per particle as a 
function of the temperature for a multi-chain configuration. 



and for different values of Vq, are sketched in Fig. 13. 

The behavior in the single chain configuration shows 
substantial differences from the multi-chain one. First 
of all, for weak values of V^', the conductivity is an in- 
creasing function of T' in the single chain case, while it 
is decreasing in the multi-chain case. This means that 
when the number of particles is small, that is the elec- 
trostatic interaction is not too strong, the thermal mo- 
tion aids the particles to overcome the potential barrier 
and does not act as a disturbance, while in the case of a 
large number of particles, the possibility of overcoming 
the barrier is sustained by the electrostatic repulsion and 
the thermal agitation is a dissipative factor. This is an 
ulterior confirmation that the dynamics of the system un- 
der investigation is a very complex interplay of driving, 
electrostatic interaction, repulsion from the constriction, 
thermal fluctuations and confinement. 

In general, in a classical model of conduction the con- 
ductivity is expected to be a decreasing function of tem- 
perature. For large values of the constriction barrier 



height either in the case of single chain and multi-chain 
structures, we found that the conductivity is not a mono- 
tonic function of temperature, although it shows a de- 
creasing trend. The presence of structure in the v x /f 
curve vs T" can be explained by the fact that with in- 
creasing temperature some of the channels formed for 
high Vq can collapse together to form new channels, as 
mentioned in previous subsection. This is confirmed by 
the analysis of the configurational energy per particle as 
a function of temperature, which is plotted in Fig. 14. 

The average potential energy exhibits the same tem- 
perature dependance as the conductivity, while in the 
case of a weak constriction barrier height the configura- 
tional energy per particle increases linearly with temper- 
ature. Another factor, which is also responsible for that 
non monotonicity, is the fact that, in the case of high bar- 
rier values, a density gradient is present and the melting 
is not homogenous, thus some parts of the system can 
be in the liquid state while others are still in the solid 
state, giving rise to complex phenomena in the transport 
properties. 
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FIG. 15: Trajectories of a system of N=560 particles for a 
high value of the driving force. The five channel structure is 
attained over the whole wire (cfr. Fig. 7(j)). 



V. OTHER DYNAMICAL PROPERTIES 

For high values of the driving force the system shows 
the phenomenon of dynamical reordering. When the driv- 
ing force is large the system, even in case of a high con- 
striction barrier, can flow in an ordered channel struc- 
ture. This is shown in Fig. 15. It is interesting to com- 
pare the trajectories of Fig. 15 with the one of Fig. 6(j). 
Above the depinning threshold (Fig. 6(j)) the channel 
structures are not homogenous, with increasing the drive 
(Fig. 15) an ordered moving structure is reached again. 

This is a well known phenomenon. Indeed, it was ob- 
served experimentally for vortex lattices in type II su- 
perconductors [5(| and for colloids [51j. The interplay 
between dynamical reordering and melting in mesoscopic 
channels was recently studied experimentally in the case 
of vortices, providing the first conclusive evidence for 
a velocity dependent melting transition |52|. The dy- 
namical reordering was also investigated theoretically for 
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FIG. 16: The ratio between potential and kinetic energy per 
particle as a function of temperature. 



CDW systems 0|. The dynamical reordering for such 
systems originates from the fact that the applied driving 
force tilts the pinning potential thereby reducing the pin- 
ning strength. When a large enough force is applied, the 
particles depin and then flow quite orderly. The same 
mechanism is responsible for the dynamical reordering 
of the studied system, with the difference that the tilted 
potential is in this case the constriciton potential. 

It is worth to study the ratio between kinetic and po- 
tential energy, averaged at every simulation step, as a 
function of the temperature. From Fig. 16 it is evident 
that the kinetic energy increases faster than the configu- 
rational energy with temperature, which is an expected 
result. What is interesting is the fact that the fitting 
curve is of the type y = c/(l + dx). The fit is excellent 
with very small errors in the fitting parameters c and d. 

Finally, we also investigated the distribution of the ve- 
locity function of the distance from the x axis, or 
in other words we studied the velocity for each conduct- 
ing channel when the system flows orderly. Motivated 
by the experimental findings of Ref . £| , that for a chain 
system of electrons on liquid helium, the particles in the 
external chains have higher velocity than the particles 
in the internal chains, we tried to clarify whether in our 
system a similar behavior occurs. Our findings are in 
contrast with the one of Ref. , actually we found that 
the internal chains have on average a velocity which is 
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5% higher than the external chains. However, this dis- 
crepancy can be explained by the circumstance that the 
pinning mechanism are different for the two systems: the 
coupling between electrons and ripplons in case of Ref. 
Q and the constriction potential in our case, and by the 
fact that the confinement potentials are not exactly the 
same in the two systems. 



VI. COMPARISON WITH OTHER DRIVEN 
SYSTEMS 

A rather general theory of periodic structures in a ran- 
dom pinning potential under the action of an external 
driving force was developed by Le Doussal and Giamarchi 
|22| . Their findings were that the periodicity in the direc- 
tion transverse to the motion leads to a different class of 
driven system: the moving glasses, with the decay of the 
translational long-range order as a power law. Similar 
considerations can be made for our system as well, but 
with the important difference that because of the confin- 
ing potential and the constriction potential the period- 
icity is broken both in the x and y-direction. For weak 
constriction barrier height long-range translational order 
is present in the x-direction, but it is softened when the 
system is moving and the temperature is increased. For 
infinite systems one of the consequences of periodicity in 
the transverse direction to the motion is that particles 
flow along static channels for uncorrelated and weak dis- 
order and that there are barriers to transverse motion. In 
our confined system the barriers to transverse motion are 
an effect resulting from confinement instead of periodic- 
ity. Most of the studies on driven lattices or glasses show 
that, at finite but low temperature, the channels broaden 
and strong non-linear effects exist in the response to the 
applied drive, though the asymptotic behavior is found 
to be linear, which is, indeed, what we found as well. 

In infinite moving systems with random pinning cen- 
ters, depending on the strength of that disorder, two 
kinds of flow are possible: (i) the elastic one, where all 
the particles move keeping their neighbors, and (ii) the 
plastic one, where part of the particles are moving in 
river-like or filamentary structures and part is pinned. 
There is a sharp crossover from the elastic to the plas- 
tic flow, related to an order-disorder transition. In our 
system we found that two kinds of flow are possible: (i) 
the elastic flow, where all the particles move orderly and 
the nearest neighbors are preserved, and (ii) the quasi- 
elastic flow, where all the particles move together, but 
creating a complex net of conducting channels, for which 
the neighbors are not kept. We also found a continuous 
crossover from the elastic to quasi-elastic flow. It is im- 
portant to stress that this difference is closely related to 
the different pinning potentials considered. To be more 
precise, in the case we investigated the pinning is due to a 
constraint rather than an actual pinning potential. That 
is the reason why we did not observe plastic flow, because 
particles cannot be strongly attracted and pinned by any 



pinning center. 

It is remarkable that for our system we also found that 
in the case of elastic depinning the velocity vs driving 
force curve scales as v oc (/ — fc) 13 , with (3 ~ 0.66, which 
is in agreement with most of the theoretical and nu- 
merical works for infinite systems exhibiting elastic flow, 
where (3 = 2/3. Thus, the value of the critical exponent 
(3 = 2/3 seems actually the signature of elastic depin- 
ning independently of the presence of confinement and 
the type of pinning potential. Naturally, in order to af- 
firm this definitely more investigations are required with 
different topologies and potentials. Furthermore, in the 
case of quasi-elastic depinning we found a critical expo- 
nent (3 ~ 0.95, which is an intermediate value between 
the case of elastic and plastic flow, where the experimen- 
tal findings give (3 ~ 2. This leads us to the conclusion 
that the quasi-elastic depinning is an intermediate regime 
between elastic and plastic depinning. 

Finally, for the elastic regime the previous theoretical 
investigations followed essential two approaches: (i) elas- 
tic theory with renormalization group techniques |22l Esj 
and (ii) perturbation theory in 1 jv [53L |5J, |55j . The first 
one explains the flowing channel structures and their mu- 
tual interactions, while the second one elucidates the v—f 
characteristics and the criticality in the depinning. De- 
spite the number of experimental and numerical data a 
detailed theoretical understanding of plastic motion re- 
mains still a challenge ■ 



VII. CONCLUSIONS 

We studied the ground state and the dynamical prop- 
erties of a classical Q1D infinite system of particles in- 
teracting through a Yukawa-type potential and with a 
Lorentzian shaped constriction potential. The system is 
confined in one direction by a parabolic potential. By MC 
simulations we found that at T = the particles arrange 
themselves in a chain-like system, where the number of 
chains are a function of the number of particles, i.e. the 
density. Depending on the height and on the interaction 
range of the constriction barrier, a density gradient in 
the chain configuration is present near the constriction. 

We studied the response of the system when an exter- 
nal driving force is applied in the not confined direction. 
We performed Langcvin molecular dynamics simulations 
with periodic boundary conditions in the not confined di- 
rection and open conditions in the confined direction for 
different values of the driving force and for different tem- 
peratures. We found that the constriction barrier and 
the friction pins the particles up to a critical value of the 
driving force. The pinned phase is a new static phase, 
with particles accumulating in the neighborhood of the 
constriction and arranging themselves in such a way to 
balance the external drive. For values of the driving force 
which are higher than the critical threshold, the particles 
can overcome the potential barrier and the system depins. 
We analyzed in detail the depinning phenomenon and 
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we found that the system can depin elastically or quasi- 
elastically depending on the strength of the constriction 
potential. The quasi-elastic flow is a new regime, where 
particles move together without keeping their neighbors. 

In the case of clastic flow the chain-like structure, 
formed at T = in the absence of external drive, is 
preserved, while in the case of quasi-elastic flow it is 
destroyed and a complex net of conducting channels is 
created. The elastic depinning is characterized by a crit- 
ical exponent, which is on average (3 ~ 0.66 and does 
not depend on the number of chains. This is in excellent 
agreement with the theoretical and numerical findings 
on 2D systems exhibiting elastic depinning. The quasi- 
clastic depinning state has a critical exponent (3 ~ 0.95. 
We demonstrated that the values of the critical exponent 
are independent of the range (i.e. screening length) of 
the interparticlc interaction. But the crossover between 
elastic and quasi-elastic flow depends on the kind of in- 
terparticle interaction. 

Furthermore, we showed that the dc conductivity is 
zero in the pinned regime, it has non-Ohmic character- 
istics after the activation of the motion and then it is 
constant, in other words the system has a non linear 
response to the applied drive. The linear regime is at- 
tained as the asymptotic behavior. The dependence of 
the conductivity with temperature and strength of the 
constriction was also investigated. We found that in the 



single chain configuration for low height of the constric- 
tion, the conductivity is an increasing function of tem- 
perature, while in the multi-chain configuration it is a 
decreasing function, as expected. For high constriction 
barrier height, the conductivity has no longer a mono- 
tonic behavior, although it has a decreasing trend. In 
these cases some structures are present in the conduc- 
tivity vs temperature curve, signaling the circumstance 
that some channels collapse or some parts of the system 
have already undergone the transition form the solid to 
the liquid state. Finally, for large values of the exter- 
nal driving force even in the case of high constriction 
barrier, the particles can flow orderly in a well defined 
channel structure, because the drive tilts the contriction 
potential, thus reducing the pinning strength, that is the 
system exhibits the phenomenon of dynamical reorder- 
ing. 
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